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Real-time Simulation of Jet Engines 
with Digital Computer (I) 

(Fabrication and Characteristics of the Simulator) 

1. Preface 

When dealing with jet engine control problems, the dynamic 
characteristics have to be determined as soon as possible. In actual 
control system tests, closed-loop tests are usually performed using a 
simulator to approximate the engine's dynamic characteristics before 
testing the actual engine prototype. This requires a simulator that 
approximates real conditions and operates in real time. 

Faster engine development times are mandatory in today's 
environment and the design and testing of control systems must take 
place during the course of development and in parallel with the 
engine development schedule. This is why high precision simulators 
have to be used in control system design so that the engine's static 
and dynamic characteristics can be simulated. Computer simulation 
uses component characteristics and engine component test data that 
are knowable in the engine design phase. 

After obtaining test results on the JRlOO lift engine, we pub- 
lished a report on methods of computing an engine's dynamic 
characteristics (shown by transfer function) from its static 
characteristics.^ The present report examines the second phase of 
those tests, tests performed in real-time simulation of engine 
characteristics and devices fabricated for that purpose. 

2. Simulating Engine Characteristics 

There are numerous reports on simulation of engine characteris- 
tics. Some simulation techniques use analog computers while others 
use digital computer methods. But the simplest use the primary delay 
shown in Figure 1 . 
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The dynamic characteristics of the engine can be expressed in the 
form: 

JNc(S)= -JWMS) 





Here 

Nc engine rpm (subscript c is the corrector) 

Wfc fuel flow 

Qc excess torque 

&2 ratio of compressor intake temperature to standard 

temperature 

^2 ratio of compressor intake pressure to standard 

pressure 

As the JRlOO example in Figure 2 shows, the gain and time con- 
stant of the transfer function are expressed by the rpm function and 
thus, they are derived as function generator rpm input and assigned 
as time constant and gain of the primary delay circuit. This simu- 
lation is extremely general, but, even then we have obtained highly 
practicable data for ground tests of gas turbines or engines 
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Figure 2 JRlOO Transfer Function 

which are components of larger systems, such as high performance con- 
trol VTOLs. In addition, all types of tests are being made both in 
Japan and abroad, and those test methods can be classified into two 
major categories, simulation by analog computer and by digital com- 
puter. 

Most representative of analog computer test simulations are 
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those conducted at the University of Michigan. Figure 3 shows a 
schematic diagram of a turboshaft engine simulation. 

The distinctive feature of this simulation is that logarithms 
are fed into the arithmetic/ logic operation circuits and it avoids 
the use of a multiplier which would be a source of noise. In addi- 
tion, devices are being constructed to closely approximate character- 
istics, such as those of compressors, using 2-input, 1-output func- 
tion generators and servos. Those devices have large configuration 
containing, overall, 60 arithmetic/ logic operational amplifiers, 10 
function generators and four 2-input function generators. Such 
large-scale devices are not always easy to use due to the need for 
adjustment, noise processing, and prevention of oscillation. Al- 
though the device is good in many respects, it has deficiencies in 
reproducibility and precision. 

The digital method is best for simulation of systems, such as 
jet engines, assembled from components, the dynamic characteristics 
of which are indicated by non-linear multi-variable functions. Digi- 
tal simulation of such systems is conducted primarily in Great 
Britain. The simulation is. used chiefly to clarify dynamic 
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characteristics, but operating time is lengthy and the simulation is 
not expressed in real time. 
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Saravanmuttoo , et. al., have simulated single- and 
dual-engine performance and examined acceleration-deceleration time 
characteristics on a compressor map locus. Because this calculation 
method repeatedly matches pressure and gas flow of an engine in 
optional state and continues its derivations in minute detail, 
computer time for single-engine simulation, reportedly, is 50 times 
real time and 200 times real time for dual-engine simulation and 
computer memory capacity is more than lOK words. 

3. Real-time Simulation of Engine Characteristics 

Requirements on control device and engine conformance are becom- 
ing stricter as engine performance increases. Satisfying the rela- 
tionship between engine input (fuel flow) and output (thrust or 
rpm's) is considered adequate with conventional, simple simulators. 
However, the conventional simulator will not meet today's stricter 
specifications. High performance and operations in real time are 
required of simulators to investigate conformity between control sys- 
tem and engine. We also have to have output on the quantitative 
status of each component. 

The following plans have been made in view of the strong demand 
for hybrid (analog/digital) simulators which operate in real time and 
can be used for control system simulation and electronic control sys- 
tem design. 




Figure 3 Simulation of Rotor Characteristics 

(1) Perforin operations in real time. 

(2) Specify simulator input such as fuel flow and environmental 
conditions in analog signals. 
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(3) Compute simulator output such as thrust, rpn's, pressure on 
all components, temperature variants and gas temperature in dig- 
ital and then output in analog. 

(4) Make simulation highly reproducible and precise. Have dy- 
namic and static characteristics of the engine conform to actual 
measurements . 

(5) Set a precision target of ±0.1% for each unit of the analog 
arithmetic/ logic operation unit used in the design of engine 
control systems and add theoretical circuits for theoretical 
design. 

(6) Use a modular mode arithmetic/logic operational amplifier in 
the analog arithmetic/ logic operation unit. Control of arith- 
metic/logic operations can also be performed by digital comput- 
er. 


(7) Use a mini-computer for overall control of the digital com- 
puter unit. 

(8) Be able to display on CRT screen, the characteristics of 
each component in the engine. 

The above is considered to be the general simulator concept. We 
plan to use a hybrid system that combines analog and digital computa- 
tion. The distribution is such that the majority of initial simula- 
tion computing is performed on a digital computer and only those 
items of integration in which the time factor enters are performed on 
analog. 

Since the actual results of setting up a program and performing 
tests showed recognizable inadequacies in simulator stability related 
to integrator precision, stability/sample and time, v;e decided to use 
a digital computer for integration items too. For I/O relations and 
for engine control system unit we used analog computer. 


6 



The simulator consists of the following units: 

(1) Digital arithmetic/ logic operation unit 

(2) Analog arithmetic/ logic operation unit 

(3) I/O interface 

(4) Display 

Figure 4.5 shows an external view of the system. 

In actually configuring the simulator, we tried, as much as pos- 
sible, to use existing components effectively. For the digital ar- 
ithmetic/logic operation unit, we used the NEAC-3200-50 which is usu- 
ally employed in high temperature turbine measuring instruments. For 
the arithmetic/ logic operational amplifier in the analog arithmet- 
ic/logic operation unit, we used a device for engine measurement and 
high performance jet engine control research. We designed and fab- 
ricated a new I/O interface in light of the goals of a real-time sim- 
ulator and because data I/O should be performed as quickly as possi- 
ble. We achieved high speed by changing over to an AD converter that 
has a 20-microsecond conversion time. This hardware will be dis- 
cussed in further detail in Section 5. 

4. Simulator Program 

Time required for computing must be known in advance in order to 
operate a digital simulator in real time. In this simulation, numer- 
ical integration is the process (process in which time is an indepen- 
dent variable) relating time. The time range from the start of com- 
puting, by data input specification, to output, is the integrating 
range of numerical integration. When the range of computing time and 
numerical integration are the same, the simulator is a real-time sim- 
ulator, in all other instances, the time range is shortened or ex- 
panded. 

The time required for digital computer calculation varies gener- 
ally according to data mode even if the same program is used. Since 
we have no elapsed timer with which to precisely know the time re- 
quired for normal computation, computer operations must be timed us- 
ing a fixed time range known in advance. Such timing can be 
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implemented by applying a pulse of fixed period to the computer in- 
terrupt line. Figure 6 shows the operations required. If interrupt 
sampling time is greater than the predicted maximum time from start 
to end of external input, calculation and external output operations, 
input will be read in a fixed time range and output obtained. 

Figure 7 shows the simulator program configuration. The two 
chief components of the simulator program are based on actual 
measured characteristics of each engine component. They are: non- 
linear simultaneous equation operations to match pressure and flow 
and first order differential equation operations which use results 
from the non-linear simultaneous equation operations to derive 
torque, perform numerical integration and calculate rpm. 
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Figure 7 Simulator Program Configuration 

The methods of solving non-linear simultaneous equations are 
numerous. All require many repetitive calculations which mean in- 
creased computer time. So, before performing simulation operations, 
we use a method of matching calculations on corrected fuel flow and 
corrected rpm combinations. We then store the resulting turbine 
pressure ratio in memory. The simulator can reference memory for 
optional input values and initial values and it can compute without 
repeatedly calculating other status variables. 


When operating the computer in real time , computer time deter- 
mines the step-size of numerical integration. When using high- 
powered computer methods (Runge-Kutta Method, Adams Method, etc.) to 
perform extremely precise numerical integration, computer time and 














computing instability greatly increase. And, if input values in the 
integration range rapidly change, abnormal values will be produced 
and when those values are used in the simulator poor conditions will 
result. Since the engine time constant, even if low, will be on the 
order of 0.4 seconds, using numerical integration (linear approxima- 
tion integration) of the oiler for a step-size of several microsec- 
onds or several tens of microseconds will give good stability and 
adequate precision. 

We will now show an example of this method where simulation was 
performed on the JRIOOH lift jet engine. Figure 8 shows the charac- 
teristics of the components in the JRIOOH that was used. We created 
two program formats, one in FORTRAN, the other in assembler language. 
The floating point method was used for arithmetic/ logic operations in 
the former, the fixed point method, in which status variables were 
normalized, was used in the latter. The ratio of computing time for 
both was 20:1. 

Figure 11 is a flow chart of simulation, and subroutines were 
created for each component category based on the data on JRIOOH com- 
ponent characteristics in Figures 8 to 10. 

The subroutines are multi-dimensional approximation expressions. 
Their values are derived by the methods described in item 4 in the 
bibliography. 
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Corrected Air Flow iCr/T./p,) 



Figure 8 JRIOOH Compressor Characteristics 



Fig. 9 JRlOOH Turbine 
Characteristics 



Figure 10 JRIOOH Nozzle 
Characteristics 


In the method of matching computation, we first postulated the 
turbine pressurization Pij/Ps, then calculated the status vari- 
ables for each component to derive turbine gas flow. Because this 
gas flow is not equivalent to the gas flow derived directly from 
P4/P5, P4/P5 was corrected to reduce the difference. This 
was repeated to derive the matching value. 
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Table 1 Component Subroutine I/O 


r, 

I 

Compressoi] c,/f, rp, 

P I 



Turbine 


Subroutine 

t 

> Input 

Output 

Compressor 

1 

; Nl^e, zt 

Gc VT/5, Ts/Tj, ric 

Turbine 

4 

! PiiPi 

74/ T5, Gif ip, 

Nozzle 

GjryPTTiPi 

yjPi 

Thrust 

' PalPi 

1 

-J 

F 


Nozzle 


p, 


p. 

. J 

F. 

rc,.T,i 

P>~ 

p> I 




p pressure 
T temperature 
G gas flow 
Gf fuel flow 
N rpm ' s 

I rotor pole moment of inertia 
Pt turbine generated torque 
Pc compressor intake torque 
h enthalpimetric 


6: ratio of standard tempera- 
ture to compressor intake 
temperature 

6: ratio of standard pressure 

to local pressure 

nc: compressor efficiency 

nt: turbine efficiency 

Gc, Gt, Gn: gas flow in each 

unit 

F: thrust 

Pa; atmospheric pressure 


Figure 11 Simulator Flow 
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Figure 12 JRlOOH Static Characteristics 

Status variables computed in this way are stored in a common 
area to which any program can refer and freely perform I/O selection 
and scaling. 

Figure 12 compares test values with static characteristics de- 
rived by the simulator. Figures 13 and 14 are engine transfer func- 
tions derived by the simulator. They are compared with previously 
announced measurement data on dynamic characteristics. Figure 15 is 
a line graph of compressor characteristics when step fuel flow is 
assigned in combination with electronic fuel control system. En- 
vironmental temperature is a parameter. Figure 16 is the data^ for 
coupling with a digitally controlled fuel control system. The solid 
lines are data from the simulator, the broken lines are data from 
actual engine use. Test conditions are somewhat different, but we 
can see a good mathc in general trend. Figure 17 shows status changes 
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in each unit of the engine when a change in step fuel flow is as- 
signed. Figure 18 is the data when the program is written in FOR- 
TRAN. Sample time is 70 microseconds. We can see the influence of 
the' sample. 
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Engine Time Constant n [SEC] 














Figure 18 Affect of Sample Time 

5. Structure and Characteristics of the Simulator 

The simulator can be divided into digital arithmetic/ logic oper- 
ation unit, analog arithmetic/ logic unit, I/O interface and display 
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unit. A schematic of the entire configuration is shown in Figure 19. 
The digital arithmetic/ logic operation unit's CPU can perform all 
simulator controls. Fixed time interval processing can be done by 
clock. Operation results are recorded by pen-write recorder, X-Y 
recorder and storage type display. A disk drive is used for tempo- 
rary data storage to store operating systems and subroutines. Each 
will be described in the following paragraphs. 

5.1 Digital arithmetic/ logic operation unit 

The digital arithmetic/ logic operation unit consists of the fol- 
lowing: 

(1) Central processing unit 

(2) Teletypewriter 

(3) Paper tape punch 

(4) Paper tape reader 

(5) Disk drive storage 

5.1.1 Central processing Unit (CPU) 

The CPU is a 16-bit control computer in parallel arithmetic for- 
mat. High speed arithmetic operation circuits have been added to 
basic commands for the fast multiplication/division required in 
real-time simulation. Four priority interrupt circuits have been 
added to the standard interrupt circuits for I/O devices. Each has 
priority over the standard interrupt circuits. One is used by the 
clock. Table 2 gives the specifications for the CPU. Figure 20 
shows the relationship between internal registers and I/O circuits. 
The A register is the bus for external output and basic operation. 

The B register is for double-length operations. The Y register is 
for specification of memory addresses. The D register is the bus for 
external input and distribution. The M register is for memory I/O. 

The X register is the index register. When an interrupt signal is 
externally input, the priority interrupt moves the program to a 
specified subroutine. It then returns to the original program after 
the second program is terminated. A priority interrupt has higher 
hardware priority than a standard interrupt. A job can be cancelled 
and a move made to a specified program even when peripherals such as 
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the teletypewriter are in operation. Both the standard interrupt and 
the priority interrupt are designed so that an interrupt inhibit and 
interrupt release can be assigned at any time by program-mask-set 
commands. Figure 21 lists the interrupt mask assignments. 

5.1.2 Teletypewriter 

The teletypewriter is used for typed output of program creation, 
measurement data, and computer results. The specifications for the 
teletypewriter are listed below 


1) 

Model 


ASR-33 

2) 

Print 

speed 

10 

cps 

3) 

Paper 

tape read speed 

10 

cps 

4) 

Paper 

tape punch speed 

10 

cps 


I Analog 
OP , 


3201 


Multi- 


er 


. Sampte 

anr 

j Hold 


AD 


AnaloE 

linB 
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OP 


( fl) 


Digital 

Output 

Circuit 


Pen 
recofcder 
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face Con- 
troller 


8CH| 
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Conv 
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]lock 


Display 


Data Signals 
Control Signals 


Figure 19 Schematic Diagram of the Simulator 
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Name 

Format 




NEAC 3200-50 


Speed 


° 16-bit parallel binary 
°Two's complementary operation 

“Current equalizer mode random access ferrite core memory, 
8K words 


Single-address mode with multi-level indirect 
and indexing 

Memory cycle time 

0 . 96psec 

Addition 

1.92 

It 

Subtraction 

1.92 

n 

Multiplication 

5.28 

n * 

Division 

10.56 

n * 

Double precision addition 

2.88 

n * 

Single-word I/O transfer 
Time multiplex 

I/O transfer 260KHz 

1.92 

at DMC* 

n 


>lMHz at DMA* 


addressing 


Power Supply 

“115V AC ± 10%, 50/60HZ ± 1.5Hz, IKW, power supply failure 
interrupt priority 
Signal Level 

“Logical 0 , OV DC 
“Logical 1, +6V DC 
Standard I/O Line 

“10-bit address pass 

“16-bit input pass 

“16-bit output pass 

“Priority interrupt 

“External control and sense line 


Weight 113Kg 

Ambient temperatures 0“C-45“C (CPU only) 



Figure 20 Computer Configuration 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

|ll 

12 

13 

14 

15 

16 


Accumulator 


Priority Interrupt (SMK *120) 


0TB Bit No. 

Device Name 

1 

Manual Interrupt 1 

2 

Manual Interrupt 2 

3 

Manual Interrupt 3 

4 

Manual Interrupt 4 


Priority Interrupt (SMK '20) 


0TB Bit No. 

Device Name 

1 

External Storage 

2 

(Digital Input 1) 

3 

Communications Controller 

4 

Indexing Attachment 

5 

(Digital Input 2) 

6 

(Digital Input 4) 

7 

(Digital Input 3) 

8 

(Disk Storage) 

9 

(Paper Tape Reader) 

10 

(Paper Tape Punch) 

11 

(Teletypewriter) 

12 

Card Reader [Punch] 

13 


14 

Line Printer 

15 

Memory Parity 

16 

Real-time Clock 


Devices in parentheses ( ) are presently in use. 


Figure 21 Peripheral Device Mask Assignment 
5.1.3 Paper tape punch 

The paper tape punch is used for high-speed input and retrieval 
of measurement data on paper tape. The specifications of the device 
are given belov;. 

1) Punch paper 8 unit paper tape 

2) Punch speed 110 cps 

3) Power supply 115V, 50Hz 


5.1.4 Paper tape reader 

The paper tape reader is used to read programs and store the 
compiler for program creation. The specifications for the device are 
given below. 

1 ) Type Optoelectronic 
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2) Read speed 

3) Power supply 


300 cps 
115V, 50Hz 


5.1.5 Disk storage 

The usual capacity of the mini-computer memory for control use 
is 4K to 8K words. When compiling source programs, the object pro- 
gram must be externally output on paper tape. Sub-routines are cre- 
ated in advance on paper tape. The only operation required then is 
storage of necessary subroutines by relocatable loader. There is 
little excess area for data storage. To overcome this disadvantage, 
a disk storage unit is set up for auxiliary memory. The disk storage 
drive used is in fixed-head mode with a capacity of 132K words. Ac- 
cess time is fast compared to movable head disk drives. Fixed-head 
drives have excellent MTBF features since there are few movable 
parts. A single interface links the disk storage drive to the CPU 
through the I/O-BUS (Figure 22). 



Figure 22 Disk Storage Device Interface 
The specifications for the disk storage drive are; 

(1) Model DATA DISC Company 1717 

(2) Capacity 132KW (1 word is 16 bits + parity) 
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(3) Disk rpm's 

(4) Average access time 


1800 

16.7ros 


Operation Unit Configuration 


o. of Comp. Used Operation Device Specs, etc. 
I i II 


Adder-integrator 









General Operational 
Amplifier 


Buffer 


Low-pass Filter 


Potentiometer 






BB1538A(I) BB3072, 
3 071 ( lU 


CEC 19-407 


BB1516 (I) , BB3305, 
064 


CEC19-301 


BB3072, 3071 (I) (II) 


CEC 19-105-2 


BB5002 


Invert 



Non-invert 

(can 

invert) 





Multiplier-divider 


Square-root Operato: 


Function Generator 


White-noise Gen. 


2 -input NAND 


4 -input NAND 


Flip-flop 


Flip-flop 



CEC 19-302 


CEC 19-303 


BB 1662 


Non-invert 
Gain 1 


Cross-fre- 
quency 5 [Hz] 


Linear 
round type 


3-line 


aervo-ser 

potentio- 

meter 


X1X2/X3 


ao >^10aiX 


CEC 19-501 (I) 
Shiquotex NE518G(II) 


BB4006 



Monostable 

Multi-vibrator 


Half-adder 


Schmitt-trigger 











































































1 


Clock Generator 


1 

^ Honeywell MC335 ,FA320 

1M,100K,10K 
IK, 100, 10, 

1 [Hz] re- 
trievable 

Chatter Shaper 

1 

5 

ST335 


1 

Pen-write Recorder 

8 

8 


Number of 
Channels 

X-Y Recorder 

2 

1 


CRT Oscilloscope 

4 

— 



(Note) I: Low-speed mode 
II: High-speed mode 


Jart 


Item 


Format 
Input Format 
Input Impedance 

Input Signal Leve[L ±10 [V] 

160 [dB] 


Open-loop Gairl 

Gain Stability 

(\ S. Temp-) 

Frequency Band Refraction 

(at 0(dbj 

Rated Output 
Output Irapedance 


BB 1538 A : BB 307.1 


BB 3072 


BB 3m 


BB 1516 


BB 3064 


Chopper ^Chopper 
tabilizetStabxli 


xze 


Chopper 
rStabxj ze 


l-sideGNI>l-sideGNDl.-sideGND 


0.5[MQ] 


0.5 [MQ] 
±10 [V] 
150 [dB] 
0.1[dB/“C]: 0.1 [dB/°C] 


15 [MHz] 

±10[V] 

±20[mA] ‘ 

5 [kl3] 


Input Voltage Drift 

(VS. TIME)| 

Input Voltag_e Drift, 

tvs. TEMP)| 

Input Conversion Noisej 
Temperature Range 

Power Supply! 


Input Voltage 0ffse|t±l5[,uV]^^^^ 


15 [MHz] 

bl0[V] 

±20[mA] 

5 [kC] 

hl0[;tVl 
(at 25[°C]) 


:l[/!V/24h] ±l[/zV/24h] 


±0.5[a 4V/°C] ^ ±0.2h«V/'’C] 


iDif f erentlDlf frntl 


-Izerl^plifier Amplif ierAmpllf ier 


6[/^V] rms 


1 [,aV] rms 


Remarks : fe§weYpSl 

Coutpug__ ; nutput 


0.5 [MQ] 
±10 [V] 
150 [dB] 
0.1 [dB/”C] 

15 [MHz] 

±10[V] 

±20 [mA] 

5 [k/?] 

±20[;tV] 

(at 25 [°C]) 

±1 [/tV/24h] 

±0.2[fiVrC] 

2 [a«V] rms 


'Dffrntl 

Input 


Dffrntl 

Input 


0.5 [MQ] ; 0.2 [MQ] 

i 

±10 [VI I ±10 [V] 

100 [dB] I 96 [dB] 

! 

0.1 [dB/T] : 0.l[dB/°C] 

I 

1.5 [MHz] i 1.0 [MHz] 


±10[V] ! 

±20 [mA] ; 

I 

5 [kC| I 

±0.5 [m\n 
(at 25[“C]) i 

±20[fiV/2ih]- : 

j 

±5[/^v/x] : 


il0[V] 

±10 ]mA] 

5[kQ] 

;0.5[mV] 

(at 25[“C]) 

:.50[jkV/ 24 h] 


:10[.«V/°C] 


4[/KV]rm> ' 10 [/zV] rms 


|(DC~10[kHz])^6[HzJ~l[kHz])j(6 [Hz]~i[kHz]j(DC~10[kHz]) ; fx:~10[kHz; 

-25~85[°C] i -25~25['’C] 
i 

I 

±15 [V] ! ±15[V] 


-25~85[°C] 
±15 [V] 

e yerf low 
xgnal 
—Output 


-25~85[‘C] i 0~60[°C] 
! 


±15[V| 


:15[V] 


bif frntl 


b 


ff rntl 
Input 

0.5 [MQ] 
±10 [V] 
86 [dB] 
0.1[dB/°C] 

5.0 [MHz] 

±10[V] 

± 10 [mA] 

5 [kQ] 

±2 [mV] 


6[/«V] rms 
-25~85[°C] 
±1,5[V] 


Table 4 Operational Amplifier Characteristics 
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5.2 Analog Arithmetic/Logic Operation Unit 

This unit is composed of two operation units, one high speed, 
the other low speed for convenience. As far as analog computers are 
concerned, both operation units are low-speed mode but one is called 
low- and the other high-speed because the former uses relays for in- 
tegrator control, and the latter uses electronic switching circuits. 

Figures 23 and 24 show external views of the analog arithmet- 
ic/logic operation unit. The following discussion is devoted to each 
component in the configuration. 

5.2.1 Adder-integrators 

The low speed mode uses relays and the high speed mode uses FET 
electronic switches to switch to and from the integrators' modes of 
reset, compute and hold. The circuits are shown in Figures 25 and 
26. This switching operation is set so that the computer can perform 
switching singly or simultaneously or so that manual switching can be 
performed. A mode for each integrator can be selected in the high 
speed mode. The circuit in Figure 27 is attached to each integrator 
to detect integrator overflow. Figure 28 is an external view of the 
integrator . 

5.2.2 Adder 

Differential input operational amplifiers are used in the add- 
er's operational amplifier. These are grounded on one side. Table 
4 gives the characteristics of those amplifiers. Circuit structure 
is shown in Figures 29 and 30. Figure 31 shows example data of a 
circle test performed on a combination of two integrators. The re- 
lation between elapsed time and error is shown in Figure 32. 

The lack of a chopper stabilizing circuit prevents modulated 
signals from being used for adder overflow circuits. Instead, a bi- 
polarity comparator was used. That circuit is shown in Figure 33. 
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5.2.3 Multiplier-divider 


The multiplier-divider uses pulse-width modulation format. The 
schematic diagram of the device is shown in Figure 34. Precision is 
a full-scale 0.1%. Figure 35 shows the frequency characteristics. 
Figure 36 gives test data. 



Fig. 24 Analog OP Unit (II) 






Figure 25 Integrator Circuit (Low-speed) 







Figure 26 Integrator Circuit (High-speed) 



\ 


Figure 27 Overflow Detector Circuit 
















Figure 30 Adder (High-speed) 


5.2 Square-root Operator 

The square-root operator has the same format as the multipli- 
er-divider. Figure 37 is the schematic diagram, Figure 38 shows the 
frequency characteristics of the operator. Test data is shown in 
Figure 39. 

5.2.5 Function Generator 

The function generator has 11 diode functions. The schematic is 
shown in Figure 40. An external view of the generator is in Figure 
41. 

5.2.6 Comparator 

Performs 2-input comparison. Output can be extracted by either 
voltage or relay contacts. Sensitivity is ImV. 

5.2.7 White-noise Generator 

Figure 42 is a schematic diagram of the white-noise generator. 

It is set up so that clock frequency can be changed to three-step and 
the upper limits of the frequency band can be changed. Table 5 gives 
the values of the power spectrum. 

5.2.8 Digital Element 

The circuit of the digital element is shown in Figure 43. The 
usable integrated circuit is DTL. The operating frequency is 2MHz 
DC. Logical 1 is +6V, logical 0 is OV. Connections can be made as 
required on the patchboard. 


33 


op-aa 




mmmmB 





mBmausmBS 









:ni 










mii 







:'. sj / ir- ?«' r.A f V 


Wctfl'ubt lr<»tf\«man*» Cep. 


Figure 31 Circle Test Data 


5.2.9 Patchboard and Control Board 


Has two surfaces one for lov.--speed use, the other for high-speed 
use. The pattern of these boards is shown in Figures 44 and 45. 


The control board is equipped v/ith-low speed and high-speed 
modes for divided use. The board is set up so that the three modes; 
single, simultaneous and computer control can be selected. In the 
high speed mode, the patch board and operational amplifiers are uni- 
fied. An external view of the device is in Figure 46. 


5.3 I/O Interface 




The I/O interface converts analog signals to digital and dig 
signals to analog between the CPU, the analog arithmetic/ logic op 
tion unit and the display. It also generates the controls signal 
for the recorder. The I/O interface consists of the following de 
vices. 
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Figure 36 Multiplier-divider Data 
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Figure 43 (a) 





















Figure 47 I/O Interface Schematic 



















Figure 48 External View of I/O Interface 
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Figure 49 Multiplexer Circuit 

(1) Multiplexer 

(2) Sample-and-hold 

(3) AD Converter 

(4) DA Converter 
(^) Controller 


The schematic diagram is in Figure 47, an external view is 
gure 48. 






5.3.1 Multiplexer 

Analog signals for fuel flow, intake air pressure, atmospheric 
pressure, etc., are the simulator's input signals. Those analog sig- 
nals are sampled by the multiplexer and input to the CPU. Modes in 
the multiplexer are changed by electronic switch and either sequen- 
tial or random-access mode can be selected. Figure 49 shows the 
switching circuits. Switching speed can be adjusted to anywhere from 
20 to 200 microseconds. Table 6 gives the specifications for the 
multiplexer. 


zero balance 

+ 15[V] -15[V] 








Figure 51 External View of AD Converter and Sample-and-hold 
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Multiplexer , / 
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Inpt Impedance 
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Temperature Range 
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racy 
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when address spec. 
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irdl'FI 
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IjgeCsampIe & 
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ri(X) [ns] 


-0.5[V]~+0.8[V] 


Temp . Range 
Power Supply 


0‘'C~70°C 
±1.5 [V] 


AD Converter (mode ADC-12Q 
(ANALOG DEVICES)) 


Convert Mode 


Output Code 


Precision 

Direct Linear-- 
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Temp Coeffi- 
cient 


/II iviJfiJ Tt 1 ■* /~rt ^ 

'•‘■’['''•'I 50 

i50 [kO] or leil. Acquire Time! cr„.,w n m ro/n 
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^ . n ^ 1 of fxnal value) 
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lonversion 

Speed 

Lnput Voltage 
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Output Level | 

Cemperature . 
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Power Supply 


Continual Comparator 

12 bit (bipolay, offset 
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lementary 
± i LSB 
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“0”<0.4[V]maxl TTL 
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+ 5(V] 



Table 6 Characteristics of the t/ 0 Interface Element 


5.3.2 Sample-and-hold 

The sample-and-hold circuit is used to precisely convert change- 
able analog signals to digital signals. We have used a circuit which 
makes the reception of alternating signals difficult by making all 
ground circuits independent for analog input and output circuits and 
for digital control circuits. The specifications of the circuit are 
described in Table 6. 


Table 7 12-bit AD Converter Characteristics 


Input Voltage (V) 

Output Voltage (V) 

0.000 

- 0.005 

0.100 

0.093 

0.200 

0.190 

0.200 

0.293 

0.500 

0.488 

0.700 

0 . G 93 

1.000 

0.991 

2.000 

1.987 

3.000 

2.983 

5.000 

4.976 

7.000 

6.973 

10.000 

9.966 


5.3.3 AD Converter 

Reducing time to the minimum is highly desirable in real-time 
simulation so that time relations, particularly time required for AD 
conversion, can be precise. We have used a modular-form 12-bit seri- 
al-comparator AD converter with 20 microsecond conversion time for 
precise conversion and ease of handling. The converter schematic is 
in Figure 50, the specifications are given in Table 6 and Figure 51 
is an external view of the converter. Table 7 lists the actual 
measured values. 

5.3.4 DA Converter 

An 8-channel DA converter has been added to the simulator so 
that simulator output can be applied to the external operator, re- 
corder and display. Input for the DA converter is binary 10-bit, 
output is a maximum of lOV, single polarity. Figure 52 is an exter- 
nal view of the DA converter. 

Table 8 shows the actual measured values for the DA converter. 
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In planning the controller we considered the possibilities of 
multi-operation of two computers for real-time simulation and gen- 
eralization, We decided on the following number of I/O channels. 

(a) Digital input 4 channels 

(b) Digital output 16 channels 


Of these, the AD converter circuit uses one digital input cir- 
cuit channel and the DA converter circuit, digital display and outp 
control circuit use 12 digital output channels. Figure 47 is a sch 
matic diagram of the circuits. Figure 53 is an external view of th 
controller. The commands in Table 9 were devised for computer con- 
trol. Corresponding hardware was designed and fabricated. 


In the following paragraphs we will discuss the AD converter 
circuits and the DA converter circuits in the controllers used by t 
simulator. 
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Figure 53 I/O Interface Controller 






Table 9 Interface Commands 


Device Name , 

Function 

Symbols 


Remarks 




Content 


' DA Converter 
and Line 
. Printer 

Outputs data if DA channel 
assignment is complete 

OTA' 041 

170041 ' 

I 

i 


Assigns DA channel 

OTA' 241 

170241 

Skips next command if. 
OTA' 041 is executed 

SKS'041, 

070041 

Skips next command if 
line printer preparations 
are complete 

SKS'141 

070141 

Directs print by line 
printer 

OCP'041 

030041 

Outputs data to line 
printer 

OCP' 141 

030141 

Digital 
Input 
No. 1 

Starts direct transfer if 
No.l flag is OFF 

OTA' 142 

170142 

HITAC-10 

and 

Digital 

Volt- 

meter 

Inputs data if No.l flag 
is ON 

INA' 1042 

131042 

1 Skips next command if No.l 
; flag is ON 

SKS' 042 

070042 

Skips next command if No.l 
flag is ON 

SKS'142 

070142 

Digital 
Input 
No. 2 

• 

• 

S Inputs data if No. 2 flag 
! is ON 

INA' 1043 

131043 

Digital 

Voltmeter 

Skips next command if No. 2 
flag is ON 

SKS '043 

070043 

Skips next command if No. 2 
flag is ON 

SKS' 143 

070143 

Assigns 200CH scanner 
channel 

OTA' 243 

170243 

Skips next command if 
200CH scanner flag is ON 

SKS '343 

070343 

Clears 200CH scanner flag ' 

OCP' 243 

030.243 

Starts digital voltmeter 

measurement * 

1 

OCP '043 

030043 

52 
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Inputs data if No. 3 flag INA'1044 131044 

is* ON 

Digital Skips next coininand if No. 3 SKS'044 070044 Y-316 

Input flag is ON 

No. 3 

Skips next command if No. 3 SKS'144 070144 

flag is ON 

Inputs data if No. 4 flag INA'1040 131044 

is ON 

, Digital Skips next command if No. 4 SKS'040 070040 

Input flag is ON 

No. 4 AD Con- 

Skips next command if No. 4 SKS'140 070140 verter 

flag is ON 

Starts AD converter after OTA' 240 170240 

assignment of 32CH multi- 
plexer 

Interrupt Places interrupt on exter- OCP'046 030046 

nal device (No.l) 

Places interrupt on exter- OCP'146 
nal device 


030146 










Figure 54 AD Converter Circuit 
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Figure 55 AD Converter Timing 
( 1 ) AD Converter Circuit 
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This circuit has functions such as multiplexer channel selec- 
tion, issuing of hold directives for the sample-and-hold circuit and 
input to the computer of generated and converted digital AD converter 
start signals. Figure 54 is the schematic diagram of this circuit, 
and Figure 54 is the circuit's timing chart. 


The multiplexer channel directives use the lower order 5 bits of 
the computer register enabling the selection of 32 channels. Of the 
AD converter's 11 output bits, computer input command enters the 
signed bit into the computer register's highest order bit and the 10 
other bits into the register's 10 lower order bits. The circuit is 
set up so that the computer can be interrupted by a conversion end 
signal from the AD converter. The computer can execute other jobs 
during conversion time. Table 10 gives an example of an AD converter 
program. 

Table 10 AD, DA Converter Program 


0001 



..EL 


oooa 



SQhr. 

ATOD. A 

0003 

00000 

0 000000 .A 

LAC 


000^1 

OOOOl 

7A 02/10 

OTA 

•2/iO 

0005 

00002 

000000 

HLT 


0006 

00003 

3/1 CO/10 


•/|0 

0007 

0000/1 

0 01 0C003 


/- 1 

000.5 

00005 

5/1 lO/iO 

I.>JA 

• lO/lO 

OOOJ 

00006 

0 01 C0005 

J.*! P 

*- 1 

0010 

00007 

lOO/iCO 

1;^'L 


0011 

00010 

100000 



0010 

00011 

0 03 00013 

i'uMA 

= 'srn 

0013 

00012 

-0 01 00000 


i\ 

001/1 

00013 

003777 




0001 

0002 

0003 

00000 

0 000000 

D 

SUBR 

REL 

DAC 

DTOAj d 

0004 

00001 

74 0241 


OTA 

'241 

OOOS 

00002 

000000 


HUT 


0006 

00003 

•0 02 00000 


LDA« 

D 

0007 

00004 

0 04 .00(^12 


STA 

JMP 

OOOS 

00005 

•0 02 'D00I2 


LDA« 

JMP 

0000 

00006 

74 0041 


OTA 

*41 

0010 

00007 

0 01 00006 


JMP 


001 1 

00010 

0 12 00000 


IRS 

0 

0012 

0001 1 

•0 01 00000 


J1P* 

0 

0013 

0014 

00012 


Ji'IP 

BSS 

END 

1 


(2) DA Converter Circuit 

The DA converter circuit uses 8 channels in the digital output 
circuit. Figure 56 is a schematic diagram of the circuit. The dis- 
play unit's X, Y-axis signals use two channels of the DA converter. 
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5.4 Display Unit 

When testing combinations of engine control systems with the 
simulator, the display unit displays engine acceleration path, etc., 
on a compressor characteristics chart. 


A storage-type CRT oscilloscope is used in the display unit and 
AD converter output is added to the X,Y axes. JRIOOH compressor cha- 
racteristics are shown, as an example, in Figure 57. The specifica- 
tions for the display unit are given here. 

(1) Model Techtronics 611 

(2) CRT 11 inch 

(3) Display 1024 bits x 1024 bits 


For display software, we created subroutines which use the same 
procedures as the software in the laboratory computer center's HITAC 
5020 X-Y plotter. Figure 58 and Table 11 show the interrelations and 
functions of these subroutines. 



Figure 56 Digital Output Circuit 






























Table 11 Display Subroutines 
ame 1 Function 


DIMENS 

Draws dimension lines 

ARROW 

Draws curved lines and places arrows at the ends 

DASHLN 

Draws dotted or solid lines between two points 

POLY 

Draws positive polygons of n angles 

RECT 

Draws rectangles 

CNTRLN 

Draws grids 

BAR 

Draws bars 

CIRCLE 

Draws circles, arcs and and helixes 

CURVEX 

Draws X polynomials 

CURVEY 

Draws y polynomials 

ELIPS 

Draws ellipses 

FIT 

Draws 3 coordinate points by parabolic 
approximation 

DASHLP 

Draws curved broken lines 

AXIS 

Draws coordinate axes 

NUMBER 

Draws floating point data in decimal 

LINE 

Draws curved lines 

SYMBOL 

Draws alphamerics, numerics and symbols 

PLOT 

Moves spot from present point to specified point 

WHERE 

Derives present point of spot 

PLOTS 

Initializes 

PLOTV 

Performs close processing 

FACTOR 

Assigns scale factor 

OFFSET 

Assigns scale factor 

SCALE 

Performs scaling 

CRTTYP 

Draws ASCII characters 

SYMBOB 

Entry only when drawing ASCII characters 

SYMBAB 

Draws ASCII characters and special characters 

SYMBOA 

Entry when drawing ASCII characters and special 
characters 

PLOTA 

Draws spot on straight line from present location 
to specified location 

SPOTON 

Spot ON 

PLOTST 

Initializes 

ERASE 

Erases display 

SPOTOF 

Spot OFF 

GRAXY 

Moves spot to specified location 


6. Conclusion 


We found excellent conformity between test data on both static 
and dynamic characteristics with the real-time simulator for engine 
characteristics described here. Simulation is very faithful compared 
to that in the analog format and operability is simple. Using assem- 
bler language we achieved a sample time of 4ms and in that range no 
^ practical problem was presented with software for real-time simula- 
tion. There is still room for improvement by speeding up the multi- 
plexer and sample-and-hold and for faster use of the interrupt 
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function. A method, as stated in Section 3, which performs matching 
computations in advance, stores the results in memory and references 
that memory for the status of an engine at any time has deficiencies 
for simulation programming. In generalizations where the number of 
independent variables is large, as with dual- and triple-engines, 
would require very large memory. At the present time, development 
has been completed on a simulation program which includes matching 
computations. An announcement cf that development is scheduled for 
the second phase of this report. 

With devices, too, the 0.1% element precision considered neces- 
sary for tests of engine control systems, has been satisfied. Fre- 
quency characteristics, along with other important factors, were ex- 
cellent. The I/O interface was designed for generalization, and it 
is possible to perform multi-processing with other computers. Using 
the simulator, we anticipate further progress in real-time simulation 
of multi-engines. 
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